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^ ■ Abstract 

An exploratory numerical study of the influence of heavy fermion doublets on 
the mass of the Higgs boson is performed in the decoupling limit of a chiral SU(2)l(8> 
' SU(2)r symmetric Yukawa model with mirror fermions. The behaviour of fermion 

and boson masses is investigated at infinite bare quartic coupling on 4^ ■ 8, 6^ ■ 12 



O 



O ■ and 8^-16 lattices. A first estimate of the upper bound on the renormalized quartic 

! coupling as a function of the renormalized Yukawa-coupling is given. 



X 



1 Introduction 

Recent LEP measurements [Q fix the number of light neutrinos to three, therefore a simple 
further repetition of fermion families is excluded. Extensions of the minimal Standard Model 
by heavy fermions are, however, possible. Examples are: a fourth fermion family with heavy 
neutrino (for a recent reference see 0), or a duplication of the three light families by heavy 
mirror families p, Some limitations on the number of heavy fermions follow from studies 
of 1-loop radiative corrections ^ because of the non-decoupling of heavy fermions. The 
question of non-decoupling in higher loop orders is, however, open. In fact, one of the goals of 
lattice studies is to investigate this in the nonperturbative regime of couplings. 

The lattice formulation of the electroweak Standard Model is difficult because of the fermion 
doublers 0. In fact, at present no completely satisfactory formulation is known 0: if one 
insists on explicit chiral gauge invariance, then mirror fermion fields have to be introduced P|, 



otherwise one has to fix the gauge as in the "Rome-approach" |jTO|. (In a recently proposed 
method |TT] a fifth extra dimension has to be introduced.) 
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There is, however, an interesting hmit of the Standard Model which can be numerically 
simulated by present techniques. Namely, if the SU(3)coiour ® U(l)hypercharge gauge couplings 
are neglected, then, as a consequence of the pseudo-reality of SU(2) representations, mirror 
fermions can be transformed to normal fermions by charge conjugation, and Yukawa models 
with an even number (Nf) of degenerate fermion doublets can be simulated [Q. (For instance, 
Nf = 4 corresponds to a heavy degenerate fermion family.) This can be done at nonzero 
Yukawa-couplings for both fermion (G^) and mirror fermion (G^), or by keeping and the 
fermion-mirror-fermion mixing mass /i^^ at zero, and thereby decoupling the mirror fermion 
exactly from the physical spectrum [IT2 . 



In the present paper we choose this second way, where the exact decoupling in the continuum 
limit is assured by the Golterman-Petcher fermion shift symmetry |13|. This symmetry is exact 
at Gv — yUvx ~ 0' implies a set of identities, which makes the parameter tuning easier. 
Another interpretation of the decoupling limit also deserves attention. Namely, interchanging 
the roles of fermion and mirror fermion by considering x to be the fermion and ip the mirror 
fermion, the decoupling scenario turns out to be a rather good approximation of the situation 
in phenomenological models with mirror fermions 0]. This is due to the fact that all known 
physical fermions have very small Yukawa-couplings. The only fermion states with strong 
Yukawa-couplings would be the members of the mirror families, if they would exist. In fact, 
the smallness of the known fermion masses on the electroweak scale could then be explained 
by the approximate validity of the Golterman-Petcher shift symmetry. 

An important set of questions for the numerical simulations is concentrated around the 
"allowed range of renormalized couplings" , which is cut-off dependent and shrinks to zero for 
infinite cut-off if the continuum limit is trivial. It is expected on the basis of 1-loop perturbation 
theory that, as a function of the renormalized Yukawa-coupling, the allowed region for the 
renormalized quartic coupling is limited by an upper bound obtained at infinite bare quartic 
coupling (A = oo), and by a lower bound called "vacuum stability bound" reached at zero bare 
quartic coupling (A = 0) (for discussions see M, [TB|, H). In the present work the numerical 



simulations are restricted to A = oo, and first results on the behaviour of the upper bound are 
obtained. The study of the A ^ limit is postponed to future work. 



2 Lattice action, decoupling, physical quantities 



Our numerical simulations were performed in the chiral SU(2)l ® SU(2)r symmetric Yukawa 
model with Nf = 2 mirror pairs of fermion doublet fields in the decoupling limit. The conven- 
tions in the lattice action and the definition of different renormalized physical quantities closely 
follow our previous papers on U(1)l ® U(1)r |15| and SU(2)l (S> SU(2)r symmetric 
models. Therefore we only repeat here the most essential formulae, and include the definitions 
specific to the present investigation. 

The lattice action is a sum of the 0(4) (= SU(2)l ® SU(2)r) symmetric pure scalar part 
and fermionic part S"^: 

S = Syy + Sqr ■ (1) 

(fix is the 2 ® 2 matrix scalar field, and \E'x = {ipx,Xx) stands for the mirror pair of fermion 
doublet fields (usually ip is the fermion doublet and x the mirror fermion doublet). In the usual 
normalization conventions for numerical simulations we have 



3^ = J2{ oTr ivtVx) + A 



1 ]2 

i^t^x) - 1 



H Tr (<^++^<^, 



2 



±4 



At=±i 

[(^iJxV'x ^Lx) + {'4^Lx¥>x'ipRx)] + [iXRxfxXLx) + iXLxVtXRx)] } • 



(2) 



Here K is the fermion hopping parameter, r the Wilson-parameter, which will be fixed to r = 1 
in the numerical simulations, and the indices L, R denote, as usual, the chiral components of 
fermion fields. In this normalization the fermion-mirror-fermion mixing mass is fi^^ = 1 — 8rK. 
In the limit A — *■ cxd the length of the scalar field is frozen to unity, therefore in only the 
term proportional to k is relevant. 

The consequence of the Golterman-Petcher identities is that at = all higher vertex 
functions containing the x-field vanish identically, and the X'X ^iid x-4' components of the 
inverse fermion propagator r^(p) are equal to the corresponding components of the free inverse 
propagator |T6|. In the broken phase the small momentum {p 0) behaviour of Txf, defines 
the renormalized ?/'-mass and wave function renormalization factor Z^, therefore 



fvp(p) = M + i-fpN + 0{p^) 



V 



(3) 



?7 ■ p 



Here /xq = /^i/>x/(2-f^) = (1 — 8rK) /{2K) and, as usual, p^ = sinp^ and p^ = 2sin|p^. The 
propagator is the inverse of F^. With the notation 



(4) 



we have 



A*(p) = f*(p)-i = A-t^-pB + Oip') = ip' + Z^iilf + lil^p'' 



Z^lfiR^p"^ - 27 ■ p(p^ + Z^fil)] Z^fipi^p"^ + Z^^l + i7 ■ p/iij^) 
Zi,fipip'^ + Z^fil + i-f ■ pfiR^) -Z^nlun^ - i-f ■ p{p^ + + Z^^l) 



+ 0(p2 



(5) 



This shows that in the broken phase near p = the elements of Aij, (p) are rapidly changing, 
unless /ip is very small. Consider, for instance, A^^: 



A 



-1 



«7 ■ p 



p2 + Z^fil + ^il^{l + Z^^ 



- + Oip' 



(6) 



One sees that in the phase with broken symmetry [fiR^ ^ 0) for fiQ and p both being small 
there is a qualitative change of the behaviour as a function of p around p^ ~ fiQ. The limits 
fio —>■ and p — *• cannot be interchanged. The correct order is to take first ^ and then 
p ^ 0. Smooth behaviour near /xq = is reached only if 



Ho = 0{p^ 



(7) 
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In this case /i^/p^ = O(p^), and the components of the propagator are 



A^^(p) = + 0(p ) , 

A,X(P) = Km = '^-P J'^'^U = 0{P) , K.{p) = --a^ + 0{p^). (8) 

P KH'Rip ~^ P ) P 

This hmit is best approximated if for a given momentum p we choose 

T 

/^P = /^o + = • (9) 

In fact, since our fermionic renormalized quantities are defined at the smallest timelike momen- 
tum Pmin (and zero spacelike momenta), we took 

(10) 

Note that this also implies A^^(pmin) = 0, which ensures that no mixing between ifj and x 
occurs. On a lattice with time extension T the smallest fermion momentum in the timelike 
direction with our antiperiodic boundary conditions is Pmin = tt/T, therefore this condition 
gives a T-dependent hopping parameter K > Kcr = l/(8r). For T — > oo one has, of course, 

K ^ K,r. 

In summary: since in the broken phase the decoupling situation is reached only for very 
small fio, in terms of the hopping parameter one has to be so close to = Kcr = l/(8r) that 
trying to perform simulations at < Kcr and to extrapolate to K = Kcr does not pay in the 
parameter region we studied (see in next section). Therefore in our simulations we always took 
/^o^/ip_=0(ir^ir,. = 0.125). 

In numerical simulations one can determine the fermion propagator A,j,. For the renormal- 
ized quantities one needs the inverse propagator Fij,. This could, in principle, be obtained for a 
given momentum by the numerical inversion of the 8 ® 8 matrix in spmoi-ip-x space, but this 



would introduce large fluctuations in the results. In [jT4|, |T5[ we used an analytic inversion up 
to O(p^). Here it is better to take 

M = (A + p'^BA-^B)-^ , 

N = A-^BM = A~^B{A + p^BA'^B)-^ , (11) 

which follows from the spinor structure A,j, = A — i'-y ■ pB, and has at most O(p^) corrections. 
In addition, for the exactly known elements of r,j, we took the values in (P), and (|ri|) was used 
only for the tp-tp component of Fij,. 

The definition of the renormalized physical quantities can be taken over in most cases 



with trivial modifications from the U(1)l ® U(l)pi, model |p!^ , |T5[. In the broken phase, where 
most of our runs were performed, we also use the "constraint correlations" obtained after an 
SU(2) rotation of the average (/^-field into a fixed ("cr") direction. (The three Goldstone boson 
components of the ip-field perpendicular to a are denoted by (a = 1, 2, 3).) In order to avoid 
infrared singularities, external vr-legs are usually set to the smallest nonzero momentum on our 
■ T lattices. 

The renormalized Yukawa-couplings can be defined in different ways. One definition is 
given by the ratio of the mass to the renormalized vacuum expectation value, as for instance 
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Gb.iI) = f^Ri(>/vR. It is interesting to compare this to the renormahzed Yukawa couphngs ob- 
tained through Goldstone-fermion-antifermion vertex functions. These renormahzed Yukawa 
couphngs, denoted as G^^^^ and G^^^^, where a = 1, 2, 3, are defined by 

. ]6>.^-,., = ^tMZ,'^^G^:Hz,'^rrniq.) , (12) 

-^75TaG^4 J V^T 

where no summation over a is apphed, and ^4, p4, are the 4th components of the momenta of 
Goldstone boson, fermion and anti-fermion, respectively. We have set the spatial components of 
all momenta to zero. The appearance of the Kronecker-delta above is due to energy-momentum 
conservation. The renormahzed 2-point fermion vertex function Tpt at small p = (0,0, 0,^4) is 
given as 



and 



f «(P4) ^^7,P. + Mn, Mn = ( ^^^^^ 

V f^R (jrxVr 



is the connected part of the iTa-"^-^ 3-point Green's function and 7ra(x) are the Goldstone fields 
{"^{y) = is the fermion field). Using the fact that the renormahzed couplings are the same 
for all three Goldstone bosons, and 



we obtain 



Troimc (75) = 4 , Tr5(/(2) (ran) = 25a 



GrI \ . _ -ikj 

_./'•^^ Ok_-n+n — ~. 7= 



-Gf^ "''-^^'^ 2AVZt 



■Y.'^^SUi2){raTTn^racb5rR{p,)Z^'^'G^:HZ^'^Y^Riqi)]} ■ (13) 
a=l 

Because of the existence of massless Goldstone bosons in the broken phase, renormahzed quan- 
tities cannot be defined at zero momentum. For instance, the connected 3-point tt-^-^^ Green's 
function has an infrared singularity on the external vr-leg. Therefore in our simulations on -T 
lattices we choose 

2tt tt tt 

= Y , P4 = , 94 = ^ • 

After carrying out all the matrix multiplications in (p^, we get Gj^], and G^^. The expressions 
are too voluminous to be displayed here. 



Numerical simulation 



We used the Hybrid Monte Carlo algorithm ||18[. This requires the fiavour duplication of the 



fermion spectrum. If the fermion matrix in the action (|^) is denoted by Q, then the replica 
fiavours have Q^, therefore for them the roles of ip and x a^i'fi interchanged: x is the "fermion" 
and ip the "mirror fermion" . Since in the model under consideration the fermions are equivalent 
to mirror fermions, in the decoupling limit G^ = = the model describes two degenerate 
fermion doublets (corresponding to the z/^-fields) and two massless "sterile" doublets (belonging 
to the x's), which have no interactions with the physical sector. 
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Table 1: Comparison of matrix inversion algorithms for various 5o and G^. 
8^ ■ 16 lattice, = 0.1, k = 0.15, K = 0.1, A = 1.0 . 



^0 




CGA 


MRA(l) 


MRA(2) 


MRA(3) 




n n 




o.u 




o.o 




n n 

u.u 




fi 1 

U. 1 


7 n 


o.o 


10-20 


0.0 


16.8 


7.9 


8.9 


10.7 


10-25 


0.0 


21.1 


9.7 


10.7 


12.6 


10-30 


0.0 


25.2 


11.8 


12.5 


14.9 


10-8 


-0.1 


5.8 


3.5 


4.1 


5.1 


10-15 


-0.1 


10.5 


5.5 


6.5 


7.7 


10-30 


-0.1 


20.5 


10.4 


11.1 


13.0 


10-8 


0.2 


11.4 


5.0 


6.1 


7.8 


10-15 


0.2 


24.4 


8.8 


10.2 


12.2 


10-30 


0.2 


52.3 


17.6 


19.0 


22.0 



The commonly used algorithm for the inversion of the fermion matrix is the conjugate 
gradient algorithm (CGA). Motivated by the paper of Gupta et al. |19| we were also testing 
matrix inversion by minimal residual algorithm (MRA) with odd-even (o-e) decomposition in 
our SU(2) symmetric Higgs- Yukawa model with mirror pairs of fermion fields. One difference 
between our case and that of [0 is that we have a nontrivial Qee and Qoo instead of their 
matrix Ml. But the requirements for a successful implementation of the MRA are satisfied 
in our case too. Namely we can invert Qee and Qoo explicitly with an amount of time that is 
negligible in comparison with the algorithmic inversion of the full matrix Q. 

To compare the MRA with the CGA we solve the equation Q^Qp = v for some scalar field 
configuration at different values of Yukawa-couplings, hopping parameter, lattice size and 
different convergence parameter. This last quantity is defined by: 

,_\Q'-Qp-v\[_ ^^^^ 



\v 



2 



Using the CGA the solution p is accepted as soon as 6 is smaller than some prescribed Sq. The 
solution by MRA is done in two steps, first solving Q^p = v and then Qp = p, both with a 
bound for 6 that is a factor of 100 smaller than 60- 



We were also testing the "polynomial preconditioning" for the MRA, described in ||19 
Preconditioning of order n, denoted by MRA(n), is characterized by the fact that the matrix, 
which is to be inverted, contains the fermion hopping parameter K to the power 2n. The tests 
were performed at A = 1.0 and A = 10-^. No important differences were observed for different 
A. As a few test runs showed, at A = 00 the algorithm behaves similarly to A = 1.0. 

Our results are summarized in tables |I| to ^ where the CPU time necessary for solving 
Q^Qp = f is given in seconds. In these tests we could not find any gain using overrelaxation. 

Table |I] shows for = 0.1 that the smaller 5q becomes, the better is the MRA compared 
to CGA. Comparing different preconditionings MRA(l) is the best. This picture changes, if 
one looks at higher K as can be seen from tables |^ and ^. Near the critical value of K simple 
preconditioning is no longer the best choice and even with an optimal preconditioning of MRA 
the CGA performs better. 
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Table 2: Comparison of matrix inversion algorithms for various K. 
43 . 8 lattice, 60 = 10-^ = 0.3, = 0.0, k = 0.09, A = lO"*^ . 



K 


CGA 


MRA(l) 


MRA(2) 


MRA(3) 


0.10 


0.59 


0.24 


0.30 


0.37 


0.11 


0.87 


0.35 


0.41 


0.59 


0.12 


1.51 


0.91 


1.14 


1.41 


0.123 


1.70 


2.41 


2.87 


2.80 


0.124 


1.74 


4.37 


4.45 


3.90 



Table 3: Comparison of different preconditionings at K = 0.125. 
43 • 8 lattice, 60 = 10"^, G^ = 0.3, G^ = 0.0, k = 0.09, A = 10-*^ . 



CGA 


MRA(n) 


n=l 


2 


3 4 5 


6 


7 


1.87 


11.30 


8.58 


6.11 


3.95 


2.29 


12.90 


15.00 



Table | shows a comparison between MRA and CGA at = 0.1 for different lattice sizes. 
One observes that the gain of using MRA increases with the lattice size. 

Our conclusion is that in the case of small K {K < 0.1) the MRA has to be preferred. 
For some choices of parameters, e.g. in the symmetric phase, the gain by using MRA can be 
so large that it would be advantageous to do calculations at different small values of K (e.g. 
0.1,0.11,0.12) and then extrapolate to the value under investigation (e.g. 0.125). However for 
K deep in the broken phase the CGA seems to be still the best choice. For the investigation of 
the decoupling limit {K ~ Kcr = 0.125) in the broken phase it was necessary to be very close 
to Kcr, and we always used CGA. 

The technical advantage of using the decoupling method is that the number of tuned param- 
eters is less, because the fermion hopping parameter is fixed ai K = Kcr- At A = 00 and at a 
fixed value of the bare Yukawa-coupling G^ one has to tune only the scalar hopping parameter 
K. The difficulty is that the presence of massless x-fermions slows down the convergence of the 
fermion matrix inversion. We also tried to improve on this by using the free fermion propagator 
for preconditioning in momentum space. In this way the number of iterations in CGA can be 
reduced by a factor not larger than two. On the other hand, the computer time required for 
performing the necessary Fourier transformations is so large that the gain is completely coun- 



Table 4: Comparison of matrix inversion algorithms on different lattice sizes 



0.0, 


K = 


0.15, K 


= 0.1, A = 


1.0 . 




■T 


CGA 


MRA(l) 


MRA(2) 


43 


■8 


0.67 


0.41 


0.46 


g3 


12 


3.84 


1.97 


2.24 


83 


16 


12.60 


6.10 


7.00 
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teracted, except for very small Yukawa-couplings. This is presumably due to the "roughness" 
of typical scalar field configurations. 

Another attempt to improve the matrix inversion was to supply the CGA with an "educated 
guess" for the start vector in the iterations. This was done by means of a hopping parameter 
expansion of Q"^ up to some order n: 

\ k=l 

where D and M are the diagonal and off-diagonal parts of Q with respect to site indices, and 
n was optimized for given parameters K and G^. For K = 0.1 and < 2.4 the speed of the 
algorithm can be increased up to a factor of three in this way. On the other hand, for K above 
its critical value this improvement is not applicable. 




4 Results 

The first step in the numerical simulations was to check the phase structure at A = oo and 
K = Kcr- On the basis of experience with several other lattice Yukawa models fl^, and 



our own previous work [^, this is expected to possess several phase transitions between the 
"ferromagnetic" (FM), "antiferromagnetic" (AFM), "paramagnetic" (PM) and "ferrimagnetic" 
(FT) phases. The resulting picture in the (G^, fi;)-plane is shown in fig. 1. Due to CPU-time 
limitations we did not try to disentangle the details of the structure near the meeting point of 
the four phases, neither did we follow the shape of the Fl-phase for very strong bare Yukawa- 
coupling beyond G^ = 1.5. 

The physical phase is FM with spontaneously broken chiral symmetry. Therefore we fixed 
Gjjj = 0.3, 0.6, 1.0 and performed a series of runs in the K-ranges shown in fig. 1 by the 
dashed lines. Most of the time 4^ ■ 8 and 6^ ■ 12 lattices were taken. In a few particularly 
important points, for instance at (k = 0.27, G^ = 0.3) and (k = 0.15, G^ = 0.6), in addition 
to 4^ ■ 8 and 6^ ■ 12 also an 8'^ ■ 16 run was performed. The typical run consisted of about 
1000-2000 equilibrating and 4000-10000 measured HMC trajectories. The length of trajectories 
was randomly changed by the number of classical dynamics steps between 3 and 10. The step 
length was chosen such that the average acceptance rate per trajectory stayed near 0.75. 

On our lattices the expectation value of the average of the scalar field in the cr-direction 
V = (cr^) = {4>Lx) (in short "magnetization") has a smooth behaviour across the physically 
important PM-FM phase transition (part of our data is shown in fig. 2). Furthermore, the 
magnetization always decreased with increasing lattice size, in the same way as in case of the 
pure 0(4)-symmetric 0^ model. This agrees with the expected second order phase transition, 
which is well suited for the definition of a continuum limit. 

In all data points the mass of the mirror fermion, /i/j^ and the renormalized mixing mass, fiR 
were consistent with zero within errors, in agreement with the consequences of the Golterman- 
Petcher relations. 

The behaviour of the fermion mass fiR^ and Higgs-boson mass in lattice units is shown in 
figs. 3 and 4, respectively. The fermion mass (/iij^) is decreasing monotonically to zero, as 
one approaches the phase transition from the FM-side (decreasing n). This agrees with the 
expectation, since in the PM phase, due to ~ Kcr, the fermion mass is nearly zero (at 
K = Kcr on an infinite lattice it would be exactly zero). At the same time the fermion mass is 
increasing with G^, in such a way that at = 1.0 within our limited computer time we were 
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unable to find points with really small masses. This could presumably be cured by investigating 
more points on larger lattices. In fact, the fermion mass shows in this point strong finite size 
effects, implying smaller masses on larger lattices. 

The masses of the doubler fermions for both ip and x at nonzero corners of the Brillouin-zone 
were also determined and turned out to be always above 1.5, with a slight decreasing tendency 
for increasing G^. 

The values of the Higgs-boson mass {mjia^ = rriL) were determined by fitting the constraint 
correlation in the a-channel by a form cosh()+ const, in the range 1 < t < T/2. The dependence 
of m/jo- on K in fig. 4 shows a gradual decrease and then a sharp increase for decreasing k. On 
larger lattices the values are smaller, but there is a substantial increase with increasing if the 
lattice size is kept fixed. This is depicted in fig. 5, where the averages of a few points with lowest 
Higgs-boson mass are shown. This figure also displays the strong finite size effects present on 
these lattices: in the limit of infinitely large volumes the minimum of m/j^ is expected to be zero 
at the second order phase transition between the FM-PM phases. A plausible interpretation 
of fig. 5 is that the finite size effects become stronger for larger G^, because the renormalized 
couplings become stronger. In any case, the large values at the minima represent a difficulty 
for the numerical simulations in the critical region, because large lattices are needed. Our 
experience shows (see also table that reasonably small masses mpi„ ~ 0.5 — 0.7 can be 
achieved at G^ = 0.3 on 6^ • 12, at G^ = 0.6 on 8^ ■ 16 lattices. Presumably at G^ = 1.0 lattices 
with spatial extension of at least 16^ are necessary. In general, for a physical interpretation of 
the results in the broken phase on a lattice of given size one should stay with k above the value 
where rriRcr takes its minimum. This expectation is strengthened by the comparison of 6^ ■ 12 
and 8'^ • 16 results at (k = 0.27, G^ = 0.3), which are both in the broken phase, and within 
statistical errors show no finite size effects of the renormalized couplings (see below). 

The behaviour of the a- (Higgs-boson) and vr- (Goldstone-boson) inverse propagators as a 
function of momentum is shown at G^ = 0.6, n = 0.15 in fig. 6. The method of measurement 
is the same as in Ref. [E^]. The Goldstone- (vr-) and a-propagators in momentum space are 
defined by 

Gn{p) = ^^^^^7r„x7r„yexp{ip- (x-y)}^ , (15) 

Ga{p) = /^-^Y^a^ayexpiip- {x -y)}^ . (16) 

In order to limit computer time and storage we actually measure G{p) for only one of multiple 
4- momenta p giving degenerate p^ according to the assumption that, at least for small momenta, 
G{p) is just a function of p^. The values are blocked during the MC runs with a block length 
of typically 100 configurations. The error bars are estimated with the jackknife method. As 
one can see, in this point the inverse vr-propagator extrapolates reasonably well to zero for zero 
momentum. The extrapolation of the inverse a-propagator to zero gives a Higgs-boson mass 
= 0.98±0.05, in good agreement with the value obtained from a fit of the time-dependence 
of timeslices. The curvature of the inverse propagators at this (k, G^) value is not strong. This 
allows a reasonably accurate determination of the renormalized quantities by the formulae in 



151 , assuming a linear dependence between zero and the lowest nonzero momentum. For larger 
values of the momentum the propagators are quite smooth, therefore the effect of heavy ip- 
fermion doublers is not strong. Closer to the phase transition the curvature of the inverse 
propagators near zero momentum becomes stronger, and the 7r-propagator starts to show an 
increasingly nonzero mass. We interpret the latter as a finite size effect. 
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Taking the Goldstone-boson field renormalization factor Z^^ = Zj- from the vr-propagator, 
one can determine the renormahzed scalar vacuum expectation value vr = {ax)/y/Z^, which in 
turn gives the renormahzed quartic- (gji) and Yukawa-couplings {Gr^) by 

„ _ _ (IR^ 

9r = — 2 — ' ^Ri^ = • v-"- ' J 

Vr 

These are shown in figs. 7 and 8. Although the errors are quite large, one can observe a strong 
increase of qr for decreasing k. In the k region above the minimum of rriR^ there is much less 
variation. The values on these "plateaus" show a moderate increase for increasing G^. The 
renormahzed Yukawa-coupling Gr^ is rather flat as a function of n, but increases definitely 
with G^. 

Considering only the points above the minimum of niRcr on the given lattice size as being in 
the broken phase, one can make a tentative first estimate of the upper bound on the renormal- 
ized quartic coupling (or Higgs-boson mass), as a function of the renormahzed Yukawa- coupling 
(or fermion mass). Such an estimate is shown by fig. 9 together with the perturbative estimates 
based on the 1-loop /^-functions. The agreement with the perturbative results at G^ = 0.3 and 
0.6 is good, although the renormahzed couplings are quite strong, i. e. close to the tree unitarity 
bound. The 8^ • 16 points have, unfortunately, larger statistical errors. At strong coupling the 
good agreement could partly be due to a fixed point in the ratio Qr/G^^^, which implies that 
this ratio is insensitive to the cut-off. 

A few measured physical quantities in selected typical points are collected in table |^. Com- 
paring the results at G^ = 0.6, k = 0.15 on 4^ • 8, 6^ ■ 12 and 8^ ■ 16 lattices with label b, c 
and d, respectively, one can see the evolution of the finite size effects. Between c and d there is 
much less change than between b and c, but point d on the 8^ ■ 16 lattice still does somewhat 
differ from the infinite volume limit. This may explain why the corresponding point in fig. 9 is 
higher than the 6^ ■ 12 points at larger k. The situation is better if one compares points C and 
D, where the deviation of the renormahzed couplings is within statistical errors. 

A qualitative relation of the masses at the strongest Yukawa-coupling (G^ = 1.0) is that 
the mass of the Higgs-boson niR^ is roughly twice as large as the fermion mass ^R^p. Because of 
finite size effects and limited statistics it cannot be decided at present whether the cr-particle 
is a two-fermion bound state or a resonance near threshold. In the latter case, due to the fast 
decay into a fermion pair, the physical Higgs-boson could become a very broad resonance. 

An interesting question is the behaviour of the renormahzed Yukawa-coupling defined 



by the 3-point vertex function in ([131) . is smaller than Gr^jj in all points. The measured 



'Ri> 

values on 6^ ■ 12 lattice are, for instance, G^j, = 1.0 ± 0.6 at point c and G'^j, = 1.8 ± 0.5 at 



point a. Therefore the ratio 



Ss^-^ (18) 

^R-ip 

is smaller than 1. On the 8^ ■ 16 lattice, within our statistics G^^ turned out to be difficult to 

measure. An exception is point D where we obtained G)j^ = 0.74 ± 0.05. This is also smaller 
by about a factor of 2 than the corresponding value of Gr^. The deviation could partly be due 

(3) 

to the nonzero momentum value where G)^^ was extracted. 

The measured values of Gr^ were always consistent with zero within small errors, in agree- 
ment with the consequences of the Golterman-Petcher relations. 
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Table 5: The main renormalized quantities and the bare magnetization (a) = {\(p\) for several 
bare couplings G^, and values near the minimum scalar mass attainable for the given lattice 
size. Points labelled by capital letters are at G^ — 0.3, whereas lower case and greek letters 
denote data obtained for — 0.6 and G^ — 1.0, respectively. 







T 




(a) 




'"'Ra 


H'Rip 


iJR 


n 

^ Ril) 


A 


43. 




0.24 


2807('15l 


307('4l 


1 23('ll 


M(2) 


48 ('3') 


\ 09('7'l 


B 


63. 


12 


0.24 


0.146(6) 


0.18(1) 


0.73(5) 


0.21(4) 


53(16) 


1.2(4) 


C 




12 


0.27 


0.303(3) 


0.309(13) 


0.80(5) 


0.39(2) 


20(4) 


1.25(5) 


D 


83. 


16 


0.27 


0.270(2) 


0.25(1) 


0.77(3) 


0.342(2) 


31(4) 


1.35(6) 


E 


63. 


12 


0.30 


0.4391(14) 


0.400(13) 


1.17(7) 


0.55(2) 


26(7) 


1.36(6) 


a 


83. 


16 


0.12 


0.118(5) 


0.136(15) 


0.63(8) 


0.61(2) 


80(50) 


4.5(3) 


b 


43. 


8 


0.15 


0.3358(16) 


0.361(8) 


1.60(6) 


1.9(3) 


59(7) 


4.9(5) 


c 


63. 


12 


0.15 


0.248(2) 


0.25(1) 


1.14(5) 


0.67(6) 


63(10) 


2.7(3) 


d 


83. 


16 


0.15 


0.218(3) 


0.217(17) 


0.86(6) 


0.54(4) 


52(9) 


2.5(3) 


c 


63. 


12 


0.18 


0.3524(18) 


0.36(2) 


1.23(8) 


0.86(8) 


36(6) 


2.4(3) 


f 


63. 


12 


0.21 


0.4390(17) 


0.41(2) 


1.34(8) 


1.11(3) 


32(5) 


2.71(13) 


a 


63. 


12 


-0.12 


0.189(2) 


0.243(14) 


1.79(15) 


0.95(9) 


180(40) 


3.9(4) 



5 Conclusions 

The important trends seen in our numerical data on 43 . 8, 63 . 12 and 83 • 16 lattices are the 
following: 

• The phase structure at (A = 00, K = Kcr) is qualitatively the same as in other lattice 
Yukawa- models with FM, PM, AFM and FI phases (fig. 1). 

• The FM-PM phase transition at A = 00 is smooth, probably of second order (fig. 2). 

• On most of our lattices there are strong finite size effects. In particular, large lattices are 
needed in order to bring the minimum of the Higgs-boson mass on a given lattice size 
down to interesting values below 1. As our simulations show, at moderate values of the 
bare Yukawa-coupling 83 . 16 might be enough, but for large values near G^p > 1.0 one 
will need at least presumably something like 163 . 32. 

• Considering only the K-values above the minimum of the Higgs-boson mass on a given 
lattice size, where in one point we also have evidence that finite size effects are not very 
strong, we obtained a first tentative estimate of the upper bound on the renormalized 
quartic coupling as a function of the renormalized Yukawa-coupling (fig. 9) . Up to renor- 
malized Yukawa- couplings at the tree unitarity limit, which is reached near G^ — 0.6, 
this agrees well with 1-loop perturbation theory, but further investigations are necessary 
in order to check finite size effects and extend the results towards larger G^. 

• At the strongest Yiikawa-coupling Gf — 1.0 the mass of the physical Higgs-boson is 
roughly equal to twice the heavy fermion mass. This could imply that the Higgs-boson 
is a very broad resonance, which decays very fast into a heavy fermion pair. 
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• For a given lattice size the renormalized Yukawa-coupling G/?^ defined by the fermion 
mass increases more or less linearly with up to = 1.0, where it becomes almost 
twice the tree unitarity bound ~ 2.5. Gj^^ defined by the 3-point vertex function is smaller 
than Gn^ on the 6^ ■ 12 and 8^ ■ 16 lattices. As discussed above, the finite size effects are 
particularly strong for G^ > 0.6, therefore large lattices are needed for confirmation of 
the values of Gr^. 

The question of the possible infiuence of heavy fermions in the Standard Model is important 
and very interesting. By numerical simulations at A = cxd one can obtain information on the 
upper limit on the Higgs-boson mass. The extension to smaller values of the bare quartic 
coupling, in particular to A ~ gives a lower bound related to vacuum stability. 

Note added: In writing this paper we received a recent preprint of Bock, Smit and Vink, 
where the same continuum "target" theory as ours has been numerically investigated in a 
staggered fermion formulation [E^]. 
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Figure captions 



Fig. 1. Phase structure of the SU(2)l ® SU(2)r symmetric Yukawa model at A = cxd in 
the {G^, -plane. The remaining bare parameters are fixed by the conditions /ip — 0, — 0. 
Open circles denote points in the PM phase, crosses represent points in the FM phase. The 
points in the AFM and FI phases are denoted by full circles and open squares, respectively. The 
dashed lines labelled R,S,T each show the range of k used for a systematic scan of renormalized 
parameters at fixed G^. The crosses along those lines denote the k. values where the minimum 
scalar mass in the broken phase is encountered. Solid lines connect the critical values for k 
estimated from the behaviour of (a)"^ on 4^ • 8. Dashed lines around the FI phase show the 
expected continuation of the critical lines. 

Fig. 2. The square of the magnetization as a function of k. Here and in the following 
figures, points at G^ — 0.3 are represented by triangles, points at G^ — 0.6 by squares and 
points at G-^ — 1.0 by circles. Open symbols denote the 4^ • 8 lattice, whereas fiUed-in symbols 
stand for points obtained on 6^ • 12. 

Fig. 3. The fermion mass iir^ plotted versus k. for G^ = 0.3 (triangles) and G^ = 0.6 
(squares) on lattices of size 4^ • 8 (open symbols) and 6^ • 12 (fiUed-in symbols). Errorbars are 
omitted when the variation is of the size of the symbols. It is seen that larger bare couplings 
G^ in general yield larger fermion masses. 

Fig. 4. The scalar mass mua- ploted versus k. The explanation of symbols corresponds 
to fig. 3. When approaching the phase transition the scalar masses increase again after going 
through a minimum. 

Fig. 5. The minimum of the scalar masses for different G^ on 4^ ■ 8 (open symbols) 6'^ ■ 12 
(fiUed-in symbols) and 8^ • 16 (square plus vertex). The G^- values of the two points at G^ = 1.0 
are slightly shifted in the plot to give a better separation. 

Fig. 6a. The inverse propagator for the massive scalar field (open squares) and the 
massless components (crosses) plotted versus the square of the lattice momentum p. The 
observed curvature is caused by the interaction with fermions. 

Fig. 6b. The inverse 7r-propagator from fig. 6a for the first few lattice momenta. The 
inverse of the renormalization constant is in principle determined by the slope of the curve 
through the origin, and is approximated by the slope of the straight line through the origin and 
the point with smallest nonzero momentum, as shown by the dotted line. 

Fig. 6c. The inverse cr-propagator from fig. 6a for the first few lattice momenta. Extrap- 
olating the curve to zero momentum gives an estimate for m^. 

Fig. 7. The renormalized Yukawa-coupling Gr^ versus k for different bare G^. The 
explanation of symbols is the same as in fig. 2. 
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Fig. 8. The renormalized quartic coupling gji as a function of k for three different bare 
values of G^. 

Fig. 9. The renormalized quartic coupling plotted versus Gj^^. Full data points are 
from runs on 6^ ■ 12 and represent the mean values obtained from the data points C,E and 
e,f in table ^ respectively. In addition, two runs on 8^ ■ 16 are shown, namely point D (open 
triangle plus vertex) and point d (open square plus vertex). The solid and dotted curves show 
the results for the upper bound on gji computed from the integration of the 1-loop /^-functions 
for a scale ratio A/mna = 3, 4, respectively. The first value for K/rriRa corresponds to rriRfj ~ 1 
in lattice units, whereas the latter is equivalent to mRcr — 0.75. (The cut-off A is defined to be 
TT in lattice units.) 
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